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Abstract 

In the last years, the need of the colour octet state in inclusive P-wave char- 
monium decay has been firmly established. However, the implications of this in 
the corresponding exclusive reactions have not been fully recognized. We argue for 
the necessity of the colour octet in P- and higher-wave quarkonium decay. Using 
a set of phenomenologically constructed baryon wavefunctions, we consider the XJ 
decay into octet and decuplet baryon-antibaryon pair. By doing so, we subject the 
wavefunctions to a test of applicability. We show that the colour singlet compo- 
nent alone is insufficient to account for experimental measurements and only by 
including the colour octet contribution can the partial theoretical decay widths be 
brought into range of the data. By the present and earlier applications of the set 
of wavefunctions, they show themselves to be reasonable model wavefunctions at 
around the scale ~ 10-20 GeV^ 



present address 
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1 Introduction 



In studying strong interactions exclusive processes, a more reliable method is one that 
bases on perturbative QCD, which calls to mind the likes of deep inelastic scattering. How- 
ever unlike in these inclusive reactions where the final states of hadrons can be summed 
over so no explicit knowledge of any of the hadrons is required, in exclusive processes 
one is usually interested specifically only in a few final hadrons at a time, as such the 
knowledge of the wavefunctions of these hadrons is essential. The more accurate picture 
of a hadron is the Fock state expansion where it is seen as a sum of states with increas- 
ing number of constituents starting from the state with only the valence quarks and/or 
antiquarks. The probability of the hadron being in any of these states is given by the 
modulus squared of the wavefunction associated with each of these states. 

\B) = V'valcnce | ^^92 ^3 ) + ^9 I ?! ^2 ^3 d) + ^1^2%' 9^) + • • • 

WiQIqI 9i--- 9m qiQi ■ ■ ■ gNQN) + ■■■ 

\M) = V'vaic„cekrg2') + V',|gM^?) + ^JgM?^) + --- 

Wiql 9i--- 9m qiqi ■ ■ ■ qNqn) + • • • 

Because the number of constituents is unlimited, one would require, in general, far too 
much information before an exclusive process involving even just a few hadrons can be 
studied. Fortunately as shown by Brodsky and Lepage in ref. in the so-called standard 
hard scattering approach (SHSA), when a high momentum transfer or high is involved 
in an exclusive process, not only is there factorization and so the process can be partially 
calculated perturbatively, but also only the lowest valence Fock state wavefunctions of 
the hadrons are needed. The higher Fock states are suppressed by the large Q^. Thus the 
large momentum transfer acts as a filter to let past only the lowest state of the hadrons. 
It has also been proved by Duncan and Mueller 0] that these are also true in quarkonium 
decay with large timelike momentum transfer. 

Most recent developments in inclusive quarkonium physics have seen, on the other 
hand, the need of the next higher Fock state in P-wave quarkonia 0, ^, the so-called 
colour octet, where the heavy QQ pair is in a colour octet rather than the usual colour 
singlet as in the valence state. This higher state is made up of the QQ plus a gluon. This 
is special to heavy quarkonium and can be understood in terms of a suppression on the 
level of the wavefunction due to angular momentum as we will explain in later sections. 
The advance in inclusive process involving quarkonia has yet to fully bring about the 
same level of understanding in the corresponding exclusive process. Part of our goal in 
this paper is therefore to show that colour octet is not only important in inclusive but 
also in exclusive process. We will achieve this by studying the P-wave xj decay into 
baryon-antibaryon. 

As mentioned above, knowledge of hadronic wavefunctions is very important. Even 
in large momentum transfer processes, valence wavefunctions are still necessary. Since 
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bound state wave equations in QCD are too hard to solve, other ways to obtain the 
solutions or the hadronic wavefunctions must be used. One way is via QCD moment sum 
rules 0, ^, 0]. However, it has been shown in ||^ that perturbative calculation of the 
proton magnetic form factor using a range of existing nucleon distribution amplitudes so- 
obtained gives results that are, on the average, at least a factor of two below experimental 
measurements in the range of up to 50 GeV^. It is therefore clear that perturbative 
contribution is not dominant in this rather low range of . The same conclusion has 
also been pointed out in p. Nucleon wavefunctions constructed from such distribution 
amplitudes are, therefore not applicable in reactions with far from asymptotic values of 
Q^, not to mention the fact that there might be ambiguities in the distribution amplitudes 
determined this way. 

In order to obtain a nucleon wavefunction that is suitable for application at as low as 
10 GeV^. The usually neglected but nevertheless ever present soft overlap between the 
final and initial nucleon wavefunction contribution to the magnetic form factor is taken 
to make up for the difference between the perturbative contribution and the experimental 



measurements. It was shown in |]I0[ that a nucleon model wavefunction could indeed be 
constructed this way, using constraints from valence quark distribution, J/?/' decay width 
into nucleon-antinucleon pair and existing data on nucleon electromagnetic form factor. 
Our motivation in this regard is to show that the above construction of a model nucleon 
wavefunction and it generalization to other octet and decuplet baryons [11| provide a rea- 
sonable description of the baryonic valence Fock state wavefunctions at around 10 GeV^. 
To do that we must apply these wavefunctions to exclusive processes. Our investigation 
will allow us to study this as well as the colour octet contribution. 

Since in xj decay into baryon-antibaryon, only nucleon-antinucleon in the final state 
has been measured, our calculations will show first that colour singlet contribution alone 
is not sufficient and by including that of the octet, the partial decay width can be brought 
into reasonable agreement with experiment. Then with the generalization of the nucleon 
wavefunction to the whole fiavour octet and decuplet multiplet, we can provide predictions 
for the widths of xj decay into other baryon pairs. 

The paper is organized as follows. First we review briefiy in Sec. ^ the angular mo- 
mentum suppression in the P-wave wavefunction and in Sec. |^ give theoretical argument 
for the inclusion of the colour octet component in exclusive decay process. In Sec. 
the set of phenomenologically constructed baryon wavefunctions will be presented and 
explained before we use the improved modified hard scattering scheme (MHSA) of Botts, 
Li and Sterman |T^, |T3| to obtain the singlet contribution in Sec. ^ and |^. The results 
will be compared to experiment. In Sec. |^, our method in obtaining a wavefunction of 
the higher colour octet state of the xj and the Feynman graphs of the decay process will 
be given. The details in constructing the octet contribution to the hard perturbative part 
Tii will be explained in Sec. |[ The actual calculation and the results are in Sec. ^ and 
p!0| , respectively. In the appendices, all essential details of the calculation, the propagators 
and the numerators of the graphs used in Sec. ^ are given. 
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2 Angular Momentum Suppression Of Charmonium 
Wavefunctions 



In this section, we show that there is a suppression on the charmonium wavefunction 
due simply to the angular momentum of the heavy quark-antiquark system. This has 
important consequences as we will see in the next section, where we briefly review the 



more detailed argument given in |T^. Hadronic decay of charmonium is through the 
annihilation of the charm with the anticharm, which is a short distance process because 
of the heavy charmonium mass. The annihilation size L is roughly given by L ~ 
the inverse of the heavy mass M. The charmonium wavefunction, which must enter the 
decay probability amplitude, is therefore needed predominantly atL~l/M~0 for a 
heavy charmonium decay. That is ips{L) ~ ips{0), the wavefunction at the spatial origin. 
This is the case for S-wave charmonium such as rjc or J/ip. For a P-wave charmonium, the 
wavefunction in fact vanishes at the origin, so instead one has to expand the wavefunction 
around the origin. What enters the probability amplitude is actually ipp{L) ~ Lip'p{L) ~ 
I/'0p(O). Going to momentum space, the probability amplitudes for a S- and P-wave 
decay include the wavefunction of the form 

S-wave: 4^s{.L) — ^ipsik) 

P-wave: ML) ~ LtP'p{0) ^Mk) , 

respectively, where k is the internal momentum of the charmonium, the average value (k) 
is of the order of a few hundreds MeV. So it becomes clear that P-wave decay is suppressed 
by 1/M in comparison to S-wave at the level of the wavefunction. This suppression is one 
of the reasons why colour octet is necessary for a P-wave inclusive decay. We will see in 
the next section that it also provides the reason for its inclusion in exclusive decays. 



3 Comparison Of The Large Scale Dependence Of 
Colour Singlet and colour octet 

The simplest scheme to calculate charmonium partial decay width that is built on the 
solid basis of perturbative QCD is the hard scattering approach of Brodsky and Lepage 
In this scheme, decay probability amplitude can be factorized into hard and soft 
part and is given by a convolution of hadronic distribution amplitudes 0's and the hard 
perturbative part Th- For the decay of xj i^^to baryon-antibaryon, or more specifically 
nucleon-antinucleon, this is 

~ fxAxA^) ® fN<l)N{x) ® fN<pN{x) ® Th{x) . (1) 

Since the distribution amplitudes are nothing but the hadronic wavefunctions with their 
internal transverse momenta being integrated over, there is still the decay constant / 
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accompanying each distribution amplitude. While the amplitude themselves are dimen- 
sionless, the decay constants /'s carry different mass dimensions depending on the original 
hadronic wavefunctions. From the fact that the partial decay width is given by 



partial 



\M\yM (2) 



and there is only one mass scale M in the process, namely the heavy charmonium mass, 
we can use power counting on Eq. (|l]) to compare, once the mass dimensions of the 
decay constants are known, the colour singlet and octet contribution to the decay width 
and hence their relative importance to a particular hadronic decay process. For the 



determination of the mass dimension of /, we refer to [H . 

Examining Eq. (Q), the decay constant of the valence Fock state of nucleon Jn and the 
colour octet charmonium decay constant f^J are both 3-particle wavefunction. Therefore 
they must be of mass dimension two. The colour singlet decay constant f^j, on the 
other hand, is a 2-particle wavefunction so it should be of dimension one. However, the 
fact that xj ^'^'s P-wave charmonia increases this to dimension two. The only remaining 
quantity in Eq. (|1]) that has a dimension is Th, which contains a hidden power of M. This 
power must make up the right dimension for A4, which must be dimension one in view 
of Eq. . So we can now collect all the dimensional quantities in the colour singlet and 
octet probability amplitude and determine their dependence on the large scale M. We 
get 

It is now clear that both the singlet and octet contribution for the P-wave xj ^i-re weighed 
by the inverse fifth power of the charmonium mass. The colour octet, although a higher 
Fock state of the charmonium, is not suppressed by the large scale of the decay process 
in relation to the valence singlet state. 

If there were no suppression at the level of the P-wave wavefunction as explained in 
the previous section, for example in the case of J/ ip decay, the dependence of the singlet 
contribution on M would have been Therefore the colour octet contribution can 

be neglected in the decay of a S-wave but not that of a P-wave charmonium . For more 



details of the argument above, one can consult |14 . 



4 Baryon Wavefunctions 
4.1 Octet Baryons 

As mentioned in Sec. |l], most model wavefunctions constructed from QCD moment sum 
rules ^ 0, |T^, |T6|, |1^ lead to rather unsatisfactory distribution amplitudes at mod- 
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erate [H, and they certainly do not describe experimental data. A different model 
wavefunction suitable for application at such low momentum transfers was constructed in 



10 1 . This construction is based on the following form of the nucleon wavefunction, 



-K2 + *^3l)l«^-^)} (5) 

which is the most general of the nucleon wavefunction with zero orbital angular momentum 



18|. The nucleon, being an isospin doublet, is also the lowest energy state of the baryons 



and is therefore reasonable for one to assume that it has zero orbital angular momentum. 
These conditions permit the presence of only one scalar function \E', which is a function of 
the light-cone momentum fractions Xj and the internal transverse momentum fc^j of the 



nucleon. In the notations of [jlO[, it is given in terms of the nucleon decay constant, or 
equivalently the wavefunction at the origin /at, the distribution amplitude 0i23 and the 
function containing the transverse momentum dependence Qn, by 



123 



[x,k±) = ^(a;i,a;2,X3;fc_Li,fc±2,fc±; 



8V3! 



and 



(IGtt 



fN{fJ'F)4>i23{x, fip) ^n{x, fc_L) . (6) 



3^2 



i=l i=l \ ' i=l 1=1 

are the usual constrained integration measures over momentum fractions and the internal 
transverse momenta. The function fi^v is conveniently taken to be of Gaussion form 



2\2 



'■N 



X1X2X3 



-exp 



I'M 



^ ^ fc± i / Xi 



i=l 



(8) 



where cat is a transverse size parameter which was fitted together with the decay constant 
/tv to experimental data by the procedure described in [jlO| to be = 0.75 GeV~^ and 
/Af(/^o) = 6.64 X 10~^ GeV^ at the reference scale /io = 1-0 GeV. 

The distribution amplitude (pus as well as the decay constant, all have to be evolved 
to the factorization scale fip of the process in question. In terms of the eigenstates of 
the evolution equation, the effect of the evolution is to change the coefficients of each 
eigenstate by a certain power of the log of the relevant scale. The nucleon distribution 
amplitude expressed in the Appell polynomial eigenbasis is 



0123 



n=l 



3 ~ 1 ~ 

1 + -^(I^Usi^) + |0?23(a;) 



(9) 

60a;iX2X3[l + 3a;i] . (10) 
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Under a change of scale both the coefficients of the expansion and the decay constant 
are scaled by the following factors of logarithm of the relevant scale p 



2/3/3o / w /A N \ "r^Z/^o 

fNilJ-p) = fNilJ-o) 



/ ln(/io/AQCD) V^^"° .;v. , ^ ./v., ^ / In WAqcd) V"^"".. . . 



where 7„ are the reduced anomalous dimensions and i3q is the first coefficient of the f3 
function. The last line in Eq. (|l^) is the expression for at the scale fip = f^o- 

As shown in [11], the form of the nucleon wave function need not be restricted to the 
SU(2) isospin doublet. On the contrary, one can extend it to the complete SU(3) flavour 
octet. The simplest way to do this, with SU(3) flavour symmetry breaking by the heavier 
strange quark mass taken into account, is to assume the complete baryon octet shares 
the same octet transverse size parameter, a^g = 0.75 GeV and the octet decay constant, 
fss = 6.64 X 10~^ GeV^. Then the flavour symmetry breaking effects are all put into the 
octet distribution amplitudes ^fg^g and they manifest themselves as uneven distribution 
of the light-cone momentum fractions amongst the valence quarks. As shown in Jll 



introducing an additional exponential dependence on the strange constituent quark mass 
m,, of the form 



"SI 

,2 ^2 



exp -^1 (12) 



Xj 



in the distribution amplitude for each strange quark with label j suffices for the purpose. 
Using several different values of m^, representative sets of the expansion coefficients 5^*, 
of the octet distribution amplitudes, (j)f23 can be obtained. The set 3 in [|T^ , obtained with 
TUs = 350 MeV, shows the most promise and will be used in the following investigations. 

The octet baryon wavefunctions with positive helicities can then be expressed as fol- 
lows. 



\Bs, + ) =i™/[dx][dV] {^fl3l/n/r-/2^)+*fl3l/f^/M 

-K12 + <\)\fiimjT-)} , (13) 

I A , + ) = J[dx] [d'k^] {M/f23 I ul^ dl^ 4^ ) - ^i, I u^J dl^ sX^ ) 

+ K32 - ^23l) luXd^sl^)} , (14) 

where stands for the quark flavour of quark j with colour and \E'f23 is the corre- 
sponding scalar functions in Eq. of the octet baryons. The A is a slightly different 
member of the flavour octet multiplet. Being an isospin singlet, it has to vanish under the 
action of SU (2) isospin and therefore has different signs between the different wavefunction 
components. 
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4.2 Decuplet Baryons 



As we mentioned in the introduction, our interest in this investigation is the xj decay into 
baryon-antibaryon pair. Unhke the octet baryon-antibaryon which can coupled only to 
spin S=l, a decuplet baryon-antibaryon pair can coupled to S=2 as well. This leads to the 
interesting potential possibility of X2 with 5*2 = 2 to decay into a coupled decuplet baryon 
pair with the same total third component of the spin. However within perturbative QCD, 
this is not possible due to quark and gluon couple via a vector coupling. This results in the 
well known helicity conservation or helicity sum rule, which forces the outgoing baryon- 
antibaryon to have zero total helicity |T9[ . Or in other words, they must be in a total spin 



one state. For the same reason, xo decay into baryon-antibaryon is forbidden. Therefore 
for the decuplet baryons, all we need are the helicity +1 decuplet baryon wavef unctions. 

Starting from the A^^, the simplest distribution amplitude which is symmetric be- 
tween the three u-quarks is the asymptotic distribution amplitude (pAsi^)- Using this as 
the starting point, one can likewise generalized to the whole decuplet baryon multiplet 
and introduce SU(3) flavour symmetry breaking in the same manner as in Sec. [4.1| by 
using an exponential rris dependence. This again yields several sets of representative ex- 
pansion coefficients -B^^°, for the decuplet baryon distribution amplitudes (f)i23ix). They 
are listed in Table |1|. The decuplet wavefunctions can be expressed in a similar fashion 
as before 

|A++, + ) =£^/[dx][dV] ^t23\<u''-<') (15) 

\B,o, + ) =^/[dx][dV] {^?2^\fllflU2l) + '^2l'!\fl^mf2l) 



+ ^f3^2"|/lVr+/2-)}- (16) 



and now the scalar functions are 



As shown in |TT], using the assumption that the nucleon and the delta have the same 
valence Fock state probabilities, the decuplet decay constant and the transverse size pa- 
rameter can take a range of values. We take as representative values fswit^o) = 0.0143 
GeV^ and ae.o = 0-80 GeV"^ 



5 Xj Decay In The Modified Hard Scattering Ap- 
proach 

Since the Xo decay into baryon-antibaryon is forbidden by angular momentum conserva- 
tion and the not-yet-confirmed ^Pi state he cannot preserve C-parity and parity simul- 



taneously in this decay mode |22] in a perturbative approach, only xi and X2 may have 
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Bi 


B2 


^3 


B, 


B, 


A 


0.000 


0.000 


0.000 


0.000 


0.000 


E* 


-0.547 


0.182 


-0.216 


-1.081 


0.062 




0.540 


-0.180 


-0.382 


1.742 


-0.413 



Table 1: The expansion coefficients of the distribution amphtudes 0^23 of the octet baryons 
considered in the xj decay. The parameters associated with this set of coefficents are 
fB.oino) = 0.0143 GeV2 and ae^o = 0.80 GeV-^ 



finite partial decay widths into baryon-antibaryon pair. The helicity amplitudes in co- 
variant form in terms of the baryon-antibaryon spinors ub{p, A) and vb{p, A) suitable for 
our consideration are, for xi 



Ai A2 A 



ub{Pi, Ai) Bi Y vb{P2, A2) 



and for X2 



-^Ia.a 



ub{pi, Ai) 62^" vb{P2, A2) e^,.(A) 



;i8) 



(19) 



X2 



where and e^,y are the polarization vector and tensor of xi and X2, respectively and 



Bj are their corresponding decay form factors. Note that Eq. (0) and Eq. ( |T9| ) are the 
only covariant form permitted for the corresponding helicity amplitudes. Using only 7^^, 
{P2 ~ Pi)^ and g^^ to form a vector and a symmetric tensor, they are the only form 
that can be constructed which still respect helicity conservation. The decay widths into 
baryon-antibaryon are therefore 



r(xi ^ BB) 



P^.,{Mb/M^,] 



167rM 



Xi 



\ E|-^Ua" 

A's 



P^.,{Mb/M^,] 



mz. 



SnM. 



XI 



and 



T{X2 BB) 



Pp,,XMb/M, 



X2/ 



167rM 



X2 



Ai A2A 



A's 



IOttM- 



X2 



(20) 



(21) 



Since in the standard or modified hard scattering scheme, baryons are treated as massless 
in comparison to the large scale M of the process, phase space must be corrected. This is 
taken care of by the phase space factor in the above equations given by pp.s. {z) = \/l — Az"^. 

Within the SHSA, one has factorization by which the soft infrared physics is contained 
in the light-cone wavefunctions and a perturbatively calculable hard scattering amplitude 



^In practice, xo has a surprisingly large upper bound on the partial width of the decay channel in 
question Q and the he may also have non-zero partial decay width. This should be attributed to non- 
perturbative soft physics or higher twist effects since mass corrections alone should not yield such large 
width assuming that the experimental width is near the upper limit. 
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Th- In the present problem, the hard scale is set by twice the charm quark mass 2mc, 
rather than by the charmonium mass M^^, because as we mentioned before, the cc pair 
annihilates at a much smaller size than that of the charmonium. This explains the ap- 
pearance of ml in Eq. (^) and Eq. (pT]). The decay amplitudes or equivalently, the decay 
form factors are expressed as a convolution of the hadron wavefunctions and the hard 
scattering amplitude Th- Based on this approach, there are already a number of work 
on charmonium decay into nucleon-antinucleon [0, ^ ^ However, we consider 
these as incomplete for the following reasons. First, they used a number of nucleon wave- 
functions, which did not describe the correct physics at the scale of the order of M^^. 
As shown in [^, none of these wavefunctions are able to describe data of the nucleon 
magnetic form factor. Second, their treatments of as are ambiguous given that the decay 
widths depend on a^, any small changes in the value of as used will change the width 
considerably. It is therefore not difficult to obtain a width that match the experimental 
decay widths. All one has to do is to choose the right value of as- This is very arbitrary 
in our opinion. A better way is to determine the scale at which as should be evaluated 
by using the virtualities of the internal exchanged gluons. However, these virtualities in 
the SHSA depend on the light-cone momentum fractions. One will encounter problems as 
the end point regions is approached when some of these gluon virtualities drop down to 
^QCD- "will become large and the perturbative part of the SHSA breaks down. Other 
treatment of this problem such as arbitrarily freezing at some values as the virtualities 
become small or using an equally arbitrary gluon mass is not well justified. Third, as we 
have already discussed due to the development first shown in p, ^ that the contributions 
from the next higher Fock state of the P-wave charmonium, where the cc pair is in a 
colour octet, are comparable to that of the lowest colour singlet contribution because of 
suppression by angular momentum. To the best of our knowledge, colour octet contribu- 
tions have not been taken into account in most exclusive reaction involving P- and higher 
wave charmonium. In the case of the decay into light pseudoscalar mesons, this has been 
worked out recently within the SHSA in and within the MHSA in ||21|| . 

To deal with the above deficiencies of the previous calculations, the phenomenologi- 
cally constructed nucleon wavefunction and its generalization to the whole of the octet 
and decuplet baryon multiplets |1T|] should provide better baryon wavefunctions at 
around 10 GeV^. The problem with the coupling can be dealt with successfully with 
the MHSA although it complicates the calculation with the additional, but necessary 
transverse momentum dependence for a self-consistent description. The above mentioned 
end point problem in which one runs into the infrared non-perturbative region within a 
perturbative scheme is cured by the introduction of radiative corrections in the form of 
Sudakov factor which spans the energy range between the lower factorization scale 
and the higher hard scale of the process in question. The Sudakov factor with transverse 
size dependence was first calculated in [|12| and was shown to cure the problem in the 
case of the pseudoscalar meson scattering. In the case of baryons, this will provide a cure 
of the above problem only if one supplements the Sudakov suppression factor with the 
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MAX prescription for deciding the infrared cutoff scale ^ . This amounts to choosing the 
largest transverse separation scale as the infrared cutoff, which is physical in the sense 
that very long wavelength gluons cannot resolve a "small" hadron which is colour singlet 
as far as the gluon is concerned. 



6 Colour Singlet Contribution 

P-wave charmonium decays dominantly through annihilation into gluons, which is a short 
distance process set by the scale of the charmonium mass M^^. The non-perturbative 
information of the bounded system must come from and be parametrized by the wave- 
function of the xj at small spatial separation, usually taken to be at the origin. For P-wave 
charmonium, the vanishing of the wavefunction at the origin forces the substitution of 
the wavefunction there by the first derivative of the radial wavefunction 



R'piO) 



3^3 J (27r)32M^,^,,^ 



f(i) 



(22) 



where \R'p{0)\ = 0.22 GeV^/^ g^^^ ^^jl-^^{k) are functions of the internal relative momen- 
tum of the cc system and are the reduced wavefunctions of the xj, that means a power of 
k has been extracted and put into the covariant spin part of the xj wavefunctions. This is 
the form of the wavefunction commonly used for xj- The related but more general form 
of the colour singlet wavefunctions of the xj for J = 1, 2 are 



(1) \ 



\ 5ab 



5 ah r d^fc 



^^^\k)S','>{p,k) \cc;k,p) 



V3J (27r)32M^, 

and the covariant spin wave function of xi and X2 expanded up to Oiji^) are 



(23) 
(24) 



2M 



XI 



Si'\p,k) = 



X2 



(25) 
(26) 



where K ■ p = The 5ab in the above equations is to ensure the cc is indeed in a colour 
singlet state. 

Since xj are even under charge conjugation, they can annihilate into two or three 
gluons at leading order 0{af) [^. The possible types of diagrams for cc annihilation into 
three light quark-antiquark pairs are shown in Fig. |l]. As discussed in [28|, a colour singlet 
quark-antiquark system has C-parity (—1)'^+'^, so P-wave spin-1 charmonia are all even 
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kQGQQfl- SSL 




(a) 



(b) 



Figure 1: Basic graphs that could contribute to xj colour singlet decay. But actually, 
only graphs of type (a) can contribute. 

C-parity states. Since strong interactions respect C-parity conservation, the intermediate 
two or three gluons must also be even under charge conjugation. Two gluons in a colour 
singlet state is automatically in an even C-parity state so Fig. |l](a) and (b) are possible 
at the two gluon stage. However in colour space, an examination of the colour structure 
of the baryon wavefunctions Eq. (|13|) and Eq. (|16]) show that exchanges between any two 
quark lines are symmetric but in Fig. ^ (b) there is an antisymmetric three-gluon coupling 
so it is eliminated. For the decay via three gluons in Fig. |I] (c), it is possible for three 
gluons to be even under C-parity through a f"''"^ coupling, but again symmetry in colour 
space between the three light quark lines forces the three gluons to couple via d"'^'^, which 
violates C-parity conservation. We are therefore left only with graphs of type (a). These 
can further be divided into four groups. Each group can be obtained from those shown in 
Fig. 1^ by permutations of the three light quark lines. Since it is usual to treat the heavy 
cc as a non-relativistic system, the heavy quarks share the energy and momentum equally 
or in other words, the distribution amplitude of the charmonium is taken to be a delta 
function which peaks at one-half. With this assumption, the hard scattering amplitude 
Tfj can be worked out. For P-wave decays, one has to keep the relative momentum K of 
the cc system and expand the hard part around K = since only terms quadratic in K 
survive the k integration. 

Taking Fig. ^ (a) as an example, with Xi, yi the momentum fractions and fe^^ 
the internal transverse momenta of the valence quarks in the baryon and antibaryon 
respectively, the virtualities of the internal lines are 





= xiyi {Ami) - {k±i + k'^if 




= Qi- (fc± 1 


+ki,r 


G2 


= (l-a;i)(l-yi) (4m2) - (fc^ 




= 02- (fe± 1 






= 0:3^3 (4m2) - (fexs + fels)^ 




= ^3 — {k± 3 


+ fcl3)J 


Q 


= (1 -xi)y3 (4m2) - ik^,-k\ 




= Q-{k±i 




Qc 


= 2[xi{l - yi) + yi(l - a;i)] ml 


+ (fc±i + fcii)^ 


= Qc+ (fe±i 





In fact. Fig. ^ (a) and (b) give identical contributions and the similar contributions from 
Fig. 1^ (c) and (d) can be obtained by interchanging x and y. One gets, after sorting and 
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SIS. 



(a) 



(b) 



(d) 



Figure 2: Graphs of type Fig. |T| (a) can be divided further into four groups. 



integrating out K, for the hard scattering amphtude 

21072 (47r)3 ml «,(ti)a.(t2)as(t3) 



9^3 {Gi + te) + te) {G^ + «e) (Q + te) (Q, - te) 



X 1/3 
+ (x < — ^ ?/) . 



:Qc-^e) 



(28) 



To improve perturbation theory via renormahzation group, the square of the renormal- 
ization scale fin for each running couphng is set to one of the virtuahtes ti,i = 1,2,3, 
which are chosen by [l^, |li 



h = max(^i, Q„ 1/bl 1/b'l) h = max(^2, Q, I/62, V^l) ^3 = max(^3, 1/&3, V&l) 

(29) 

with ri/ = 4 and Aqcd = 220 MeV. 

These must be convoluted with the wavefunctions to obtain the colour singlet decay 
form factors. In transverse separation space and in terms of the expressions. 



Ci{x,y) = {-ly-^ y3 xi{xi - yi) 



(30) 
(31) 



the hard scattering part is 

fj^(x,?/, 6,6') = ^-^^^ as{ti)as{t2)as{t3) S'^ {bi - b[ - + b'^ 



x|^/f«(V5^(2m,)|63|) jHl^'\^il-xM2m,)\b,-b[\) 



H^"{V^i{2m,)\b[\) .J , , Ci{x,y) 

i-i l^j y) ~\ ; 

xi + yi 



(xi + ?/ + 1)(1 -xi- yi) 



{l-xi- yi){2 -xi- yi) 
13 



{C({x,y) 



Ci{x,y) 
2-xi-yi 



4 K^^\^xi{l - yi) + yi{l - xi)mM) 
71" {xi + yi){2 - xi - yi) 

{Ci(x,y) + —^^^ i) 

^ [xi + yi){2-xi-yi)^ 

2mM Cijx, y) k'^i\^xi{1 - yi) + - xi)m,\b\\) \ 
- yi) + yi{l - xi) {2-xi- yi){xi + yi) J 
+ ((x, b) ^ [y, b')) . (32) 

With this, the decay form factor can be expressed as 

exp[-5(x,y,b,6',2m,)] \\^''{x,b)^''{y,b')\\ (33) 

where aj = 1/V2, 1 for J = 1, 2 respectively, and Sudakov correction factor evaluated at 
the scale of 2mc is included in the convolution. As mentioned earlier, the presence of this 
radiative correction in the intermediate scale range, together with the MAX prescription 
for the infrared cut-off in the Sudakov factor, b = max(6i, 62, &3, &2, ^s), renders the 
whole approach self-consistent. Actually, since there are two hadrons in the final state, 
one can have a separate infrared scale for each hadron, for example b = max(6i, 62, ^3) and 
b' = max(6']^, 62, 63). However, numerically this would make no difference so we merge the 
two into one scale. The factorization scale fip in the wavefunction is then set to fip = 1/b 
as usual. 

The spin structure of the contribution from the diagrams of Fig. ^ requires the of 
the charmonium to be equal to 5*^ of the quark line that is not attached to other light 
quark lines via gluon, i.e. the first quark line from the left of the three vertical lines 
in any of the figures in Fig. |^. This permits two possibilities for the helicities, that is 
(-|-,-|-,-) and (-|-,-,-|-) of the u and d-quark lines in Fig. g from left to right. The hard 
scattering amplitudes are, however, identical for the two arrangements of helicities. This 
gives the following sums of products of the Fourier transform of the scalar functions 
"^123 of baryon-antibaryon, given in Eq. (||) and Eq. (pT]), represented in Eq. (|33| ) by 
\\<i/^ix,b)^^iy,b')\\. 



\^^^{x, 6)^^«(y, 6')|| = 2 I ^f|3(x, b)^f,%{y, b') + ^g\{x, b)^g\{y, b') 

+ (^f|3(x, b) + ^f|i(x, b))(^f|3(y, b') + ^g\iy, b') 
+ (2< — >3)| (34) 

\\i>\xM^\y,b')\\ =6 ( ^^,{xM^23{y.b') + ^,Ax.b)^2l{y.b') 



14 



J 


rW(xj^pp) (eV) 


PDG (eV) H 


BES (eV) [||] 


1 


2.53 


75.68 


37.84 


2 


16.58 


200.00 


118.00 



Table 2: Clearly, the colour singlet contributions are insufficient in explaining the exper- 
imental data of Xj decay into pp. 



+ (2 < — > 3) 



b')\\ = 3 I ^fiiix, b)^fii{y, b') + ^ii^ix, b)^g'ny, b') 
+iffiS (x, b)^fiS (y, b') + (x, 6)^|\° (y, b' 



(35) 



(36) 



The results for the colour singlet contributions in Xi X2 decay into nucleon- 
antinucleon are shown in Table ^ Q together with the experimental measurements. It 
is clear that the singlet contribution is insufficient even with the smaller results of the 
BES collaboration |^ to account for the experimental measurements. For uncertainties 
in the theoretical estimate of the colour singlet contribution, the consistency of this calcu- 
lation within MHSA, the choice of proton wavefunction and the value of the proton decay 
constant used etc., we refer to |]14[. We will discuss more about the differences between 
the experimental results in later section. The very important colour octet contribution 
will be investigated in the following section. 



7 Colour Octet Contribution 

In colour octet cc decays into baryon-antibaryon, there is a constituent gluon in the initial 
state so the C-parity arguments given in Sec. ^ no longer hold. The diagrams in Fig. |T| 
then form the bases of three different contributing groups. The other possibilities are from 
the graphs where the cc pair annihilates into a single gluon which would not be possible, 
if the pair were in a colour singlet state. These additional groups are shown in Fig. ^ 

In order to form colour singlet baryon-antibaryon in the final states, this net colour 
from the constituent gluon must be neutralized. One could allow it to enter directly into 
one of the final baryons as a constituent of a higher Fock state of the latters. However, 
within the hard scattering approach, at least one hard gluon must be exchanged between 
all the constituent partons, therefore such contribution involving the next higher Fock 
state of the baryons will be suppressed by the hard scale of the process and also by 

■'Our present numbers supersede those previously reported in pG]. 
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Figure 3: In addition to the graphs of type Fig. |T], these form the bases of further contri- 
butions in the colour octet decay channeL 
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the smaller probability of the next Fock state. In any case, bringing in a higher state 
unnecessarily will introduce additional unknown wavef unctions. It is therefore best to 
avoid it. The introduction of the colour octet in the charmonium system is, on the 
contrary, well founded and is indeed necessary as shown in p, || and as we argued earlier. 
The alternative for colour neutralization is done by attaching the constituent gluon onto 
all possible places in the diagrams of Fig. Fig. ^ and Fig. This will generate about 
9-11 diagrams for each group and the light quark lines will also have to be permutated 
and so altogether there are over two hundred diagrams. Fortunately, the tedious algebra 
can be handled completely by the computer program FORM, and one can arrange the 
constituent gluon to be attached automatically to all possible places in each basic graph 
and worked out the colour factors as well as the perturbative hard scattering part entirely 
by using this program. In fact, none of the algebra needs to be worked out by hand. 
The colour octet wavefunction of the charmonium is 



\X?!,P) = /[d^] 0S(^i,^2,^3)^S(p) \ccg;p) . (37) 



As explained in ||TT|, the distribution amplitude of the octet charmonium is conve- 



niently taken to be delta functions that peak at the light-cone momentum = z = 0.15 
for the constituent gluon and Zi = Z2 = {1 — z)/2 for the heavy quarks while the octet 
decay constant f^J are obtained from fits in pO|, |21|]. The covariant spin wavefunctions 



for the octet states are given in |2T 



S?J{P) = ^(^ + ^xJ e,.p.P^6-7'^ 

Si'J{p) = ^{^ + M^Je^.r . (38) 

As we will take the spin of the charmonium to be pointing upward without loss of gener- 
ality, we found it more straight forward to write these as 

S[T\p) = li^ + M,MPef^ -0^) (39) 
sT\p) = \{i^ + M,.,){^^^^ef^ +0e^;y) (40) 

instead. The ep is the polarization vector of the constituent gluon and e^c is the spin S=l 
vector of the colour octet cc system. Then the numerator of any contributing graphs to 
Xj decay can all be expressed in the form of 

= A+{-iyM (41) 

and the difference in the numerator between the xi ^"^^ X2 system comes entirely from 
the sign. In the appendix where we list all the contributions from the graphs of each 
group, the numerators are all in the above format. 
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To deal with the colour octet contributions, bearing in mind that the advantages of 
the dynamical setting of the renormalization scales and the built-in Sudakov suppression 
of the end-point problematic regions of the distribution amplitudes, one could use again 
the modified hard scattering scheme |T^, as we did in the singlet contribution in the 
previous section. However, these advantages are obtained at the expense of complicating 
the expressions and the calculations of the perturbative hard part Th{x, k±) by including 
the internal transverse momenta in the propagators. These transverse momenta will have 
to be integrated out subsequently by convoluting with the hadronic wavef unctions. There- 
fore the number of integration variables can be quite high especially when we are dealing 
with baryon and antibaryon which contain three constituent quarks and antiquarks even 
at the valence level and also there are many diagrams to consider. To keep things simple, 
it is advisable to return to the standard scheme so as to deal only with distribution 
amplitudes and Th{x) without the internal transverse momenta and not the more com- 
plicated wavefunctions and TH{x,k±). This is what we will do below. We will discuss 
getting the graphs and the expressions for the individual colour octet contribution to the 
amplitude in the next section within the standard scheme. 



8 Getting The Graphs And Calculating Th Of The 
Colour Octet Contributions 

In the appendices, we give the graphs and expressions of our calculation for the colour 
octet contribution. The diagrams can be divided into ten basic groups. Each group is 
based on one basic graph in which the constituent gluon from the colour octet component 
of the charmonium has not yet been drawn or inserted. Because of the restriction of 
C-parity, in the absence of the constituent gluon, these groups individually may or may 
not exist in the colour singlet contribution to the decay. In the latter case, they survive 
solely because of the presence of the constituent gluon. These graphs are those already 
presented in Fig. |l] and Fig. ||. The group associated with each basic graph is generated 
by attaching the constituent gluon to all possible places on the basic graph except on the 
initial c-quark or c-quark when they have just emerged from the charmonium. 

The numerator of the hard part Th of each diagram from each group will be given in 
the appendices, while the denominator which is essentially a product of the propagators 
in each diagram, will not be written out explicitly individually but can be derived from 
those of the basic graphs (these will be presented with each group) by following some 
simple rules that followed from the momentum flow through the graph. Given that the 
momenta of the xj? and those of the outgoing baryon and antibaryon are P^j, Pb and 
Pg respectively, the outgoing momenta of the constituents of the baryon and antibaryon 
will be assigned the momenta Pg. = x^Pb, Pq^ = ViPs with i = 1,2,3 . Here Xi and Ui 
are the momentum fraction of the i light quark line of the outgoing baryon-antibaryon 
subjected to Y^i^i = Y^iVi = ^- In the basic graphs, the constituent gluon is not included 
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so the momentum of xj shared by the charm- ant icharm is Pc = ^i-^xj ^^"^ Pc = ^'^Pxj 
with Zi + Z2 = 1. Obviously energy- momentum conservation requires P^j = Pb + Pg- 
The product of the gluon propagators in each of the basic graphs is expressed with these 
momentum configurations with 2Pb ■ Pb = (2mc)^ ^ P| ~ 0. On the insertion of the 
constituent gluon with momentum fraction z, the condition Zi + Z2 = 1 is now replaced 
hy z + zi + Z2 = 1. Therefore the heavy quark pair will have only (1 — z)P^j instead 
of the full P^j to annihilate into gluons. The consequence is that in the new graphs the 
momenta flowing along the path or lines connecting the initial heavy quark line of the 
charmonium to the insertion point of the constituent gluon must be reduced by zP^j. 
The momenta in any intermediate heavy quark line(s) which only exist in those graphs 
where the cc annihilates into two or more gluons, may also need to be so-adjusted but that 
depends on the insertion point. This momentum shift in the basic graphs brings about a 
similar shift in the momenta carried by the propagators along the same connecting path. 
In general, this shift can be done by replacing one of the three {xi, Ui) pairs in the affected 
propagators including the intermediate heavy quark one by (xj — z,yi — z). 

The momentum fractions and the quark and gluon lines are labelled as follows for 
simpler description and they apply to all basic graphs. The top (bottom) heavy quark 
line on the left is the initial charm (anticharm) quark line. The three vertical light quark 
lines running from the bottom to the top are light quark line 1 to 3 from left to right. 
Light quarks (antiquark) of the baryon (antibaryon) emerge from the top (bottom). The 
fraction Xi and Ui label that of the quark and antiquark of the i light quark line. The 
gluon line coming directly or indirectly from the charmonium and entering the i light 
quark line is gluon line i, Gi. Whereas the labelling of the vertical light quark line is 
for the entire line, each of these gluon lines connects only two vertices and ends at these 
vertices. There will be only one gluon line entering each light quark line even though each 
light quark line may be attached to more than one gluon line. In this case, only one such 
line enters the quark line and the other is leaving it. The gluon line is then labelled with 
Gi if it is entering light quark line i and not leaving it as we mentioned above. Fig. ^ 
illustrates our labelling scheme with some examples. There are however gluon lines that 
enter a 3- or 4-gluon vertex instead of a light quark line, see Group 4 and 6 for example. 
These are not labelled. The labelling is for the purpose of identifying the graphs created 
from the basic one by the insertion of the constituent gluon. The insertion point (LP.) 
of the constituent gluon will uniquely identify each so-created graph and these graphs 
will form a group based on the basic one. Taking the xi ^"^^ X2 to be spinning up. The 
allowed helicities of the three light quark lines (Ai,A2,A3) are (+,+,—), (+,—,+) and 
(— , +, +). These three helicities configurations of each graph do not necessarily have the 
same or permutation related numerator as in the colour singlet case. The presence of the 
constituent gluon makes this not possible in general. The corresponding numerators are 
listed in the tables in the appendices. 

As mentioned above, each graph created from a basic graph by adding the constituent 
gluon is identified by the insertion point where the gluon from the colour octet state of 
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^1 X2 



^1 X2 X3 





yi y2 ys 

(a) 



yi y2 ys 

(b) 



Figure 4: Our labelling scheme as applied to (a) Group 2 and (b) Group 4. 



the charmonium joins the basic graph. The group of a basic graph is created by attaching 
the colour octet gluon to all possible places on the basic graph with the exception of the 
initial heavy quark or antiquark line. Also the creation of a graph with an intermediate 
state with one gluon and nothing else is forbidden by C-parity. More explicitly it means 
that when the basic graph contains the cc annihilation into one single gluon, the graph 
with the constituent gluon attached to this gluon is not possible. Note that this kind of 
basic graphs by themselves cannot exist without the constituent gluon because of C-parity 
and so they do not contribute to the colour singlet contribution. 

The complete graphs created from the basic ones will be labelled as follows. If the 
constituent gluon attaches to one of the light quark lines, it can be with one of the three 
quarks (antiquarks) that forms the baryon (antibaryon). These graphs are labelled as Ui 
or Li indicating that the constituent gluon joins the upper (top) or lower (bottom) part 
of the quark line i. If the quark line has one or more intermediate off-shell middle seg- 
ment, which is or are separated from the upper and lower part by quark-gluon- antiquark 
vertex or vertices, to which the octet gluon attaches itself, these are labelled Mi for one 
intermediate off-shell quark line on quark line i, for example the quark line 2 of Group 1 
(sees Fig. |^), and UMi or LMi denoting upper or lower part of the middle off-shell line if 
two intermediate middle segments exist on the light quark line i, for example the quark 
line 1 of Group 8 (see Fig. |1^). Evidently, the constituent gluon can be attached to the 
gluon lines as well. These graphs are labelled with Gi (not to be confused with the gluon 
line labelling) if the new graph is created by joining the gluon to the gluon line i by a 
3-gluon vertex. If the gluon is attached instead to an unlabelled gluon line, the graph will 
be called GR as the constituent gluon attaches to the remaining unlabelled gluon line. 
Only one such graph per group is possible at maximum and these always vanish due to 
colour. We will not discuss these graphs further. Whenever a 3-gluon vertex exists in a 
basic graph, there is the possibility of attaching the constituent gluon to this vertex to 
create a complete graph with a 4-gluon vertex. These will be labelled as 4G for one such 
possibility or 4G1 and 4G2 when two such possibilities exist as in Group 9 (see Fig. 0). 
There remains attaching the gluon to the heavy intermediate heavy quark line(s) when 
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this exists. We denote this simply by graph Q for only one heavy quark line or UQ and 
LQ in the presence of two heavy quark lines, one situated above the other on the basic 
graph, meaning the gluon is attached to the upper or lower heavy quark line. This latter 
is possible only for Group 5. 

The rules for getting the product of propagators for each new graph from that of the 
basic graphs are as follows. In determining which of the pair (xj, t/i) in the propagators 
to replace by (xj — z,yi — z), it is the i pair for graphs labelled \Ji, Li, Gi, Mi, UMi 
or LMi. The replacement should only be made to propagators along the line connecting 
the heavy quark line to the insertion point whether they contain explicitly or 

not. Modification must be made to gluon or light quark propagators through the relation 
Xi = 1 — Xj — Xk and Ui = 1 — Uj — tjk where j ^ k ^ i and also to any charm propagator 
only if it contains the pair {xi,yi) explicitly. So applying this to graph Ul of Group 5 
illustrated in Fig. |(a), only the propagators associated with the gluon line Gl and the 



upper heavy quark line in Fig. 12 need this replacement. The product becomes 



{{zi+z-xi){zi+z-yi)-l/A}{2mc)^ {{z2-x:j){z2~y3)-l/4:}{2mc)^ {xi-z){yi-z){2mcy+p^ X2y2{2mc)^ xsyz{2mcy ' 

(42) 

If it is graph L2 of Group 5 in Fig. |^(b), only the propagator of gluon hne G2 needs to 
be changed and we have 

1 1 1 1 1 

{(^i-xi)(2i-?/i)-l/4}(2mc)2 {{z2-x-j.){z2-yj,)-l/ i}{2mc)'^ xiyi{2mc)^ {x2-z){y2-z){2mc)2+p^ xzyi{2mc)'^ • 

(43) 

As a third example, applying this to Fig. ^(c) that is graph L3 of Group 1 gives 

1 1 1 1 1 

(l-z)2(2mc)2 (x3-2)(l-{/i-2)(2mc)2+p2 xiyi(2mc)2 (l-x^-z){l-yi~z){2m.c)'^+p2 {xz-z){yi-z){2m^Y+p^ • 

(44) 

It must be mentioned that in the appendices where details of each group are given, the 
propagators are written down without the usual ie. This is to save on typing but it should 
be understood implicitly that a ie should be present in each propagator. This also applies 
to the expressions here. Also, we have inserted a term in the denominators of those 
propagators which carry two poles which cannot be handled by the ie prescription. The 
is understood to be the mean-squared internal momentum of the baryons = (fcj.^). It 
is used here to prevent two possible poles in any one propagator to occur simultaneously. 
This problem has been treated in the same manner in whenever the ie prescription 
failed. The actual values of depend on the baryon wavefunctions. We used p(8) = 
415.0 MeV and p(io) = 389.0 MeV for the octet and decuplet baryons, respectively. These 
values are obtained from the respective wavefunctions. 

For 4G graphs, if the unlabelled gluon line attaches to the 3-gluon vertex connecting 
gluon line Gi and Gj, then either replace {xi,y.i) with {xi — z,yi — z) or {xj,yj) with 
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(a) 



(b) 



Figure 5: Examples of complete colour octet graphs, (a) graph Ul and (b) graph L2 of 
Group 5, and (c) graph L3 of Group 1. 




(a) (b) 

Figure 6: Examples of complete colour octet graphs with 4-gluon vertex, (a) graph 4G of 
Group 4 and (b) graph 4G of Group 10. 

{xj — z, Uj — z) in the gluon propagator and do the same in the heavy quark propagator but 
in this case only when either pair of momentum fractions appears explicitly is sufficient. 
Applying this to Group 4, we get 

1 1 1 1 1 

{(22-2:2)(z2-S/2)-l/4}(2mc)^ {l-X2-z)(l-y2-z){1mc)'^+p'^ xiyi(2mc)^ X2y2{'2mcY xzyi{2mc)'^ ' 

(45) 

and to Group 10, the product of propagator become 

1 1 1 1 1 

(l-z)2(2mc)2 (l-2){l-j/i-z){2mc)2 {l~x^-z){l-yi-z){2m^Y+p^ X2y2{'2m^Y xzyz{2mc)'^ ' 

(46) 

Eq. ( ^Ij ) and Eq. (^) are fine examples of the necessity of having to make the shift 
even if either pair of momentum fractions does not appear explicitly in the light quark 
and gluon propagators. In these example, the replacement has to be made even for 
xi + X2 + = Ul + y2 + = 1 to 1 — z because the momentum fractions are completely 
hidden. For the remaining Q, UQ and LQ graphs, it is only necessary to replace zi by 
2:1 + 2; and Z2 must be left alone. 

One must not forget that the insertion of the constituent gluon will introduce an 
additional propagator (X-Prop.) to the basic graph. This is, however, not true for the 
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4G graphs but these have a factor (2mc)^ less than the other graphs and so a factor of 
l/(2mc)^ must be multiphed to them so that a global prefactor can be written down for 
all graphs. The extra propagator or the extra factor l/(2mc)^ is given in the last table of 
numerators of each group under the column X-Prop. Similarly, the structure of group 6 
is such that it is a factor (2mc)^ less than the other group from the start because these 
graphs have one less propagator. In order to keep the same overall factor for all groups, 
l/(2mc)^ has been multiplied to the product of the basic propagators of this group. This 
extra factor of course does not correspond to any propagator but is merely a compensation 
factor. It must be mentioned that not all graphs generated from the basic graphs exist. 
Some vanish for one reason or another. For these, the additional propagators are not 
given in the tables but have the entry n. g. (not given) instead. 

There exist several more groups that we have not drawn out explicitly in the appen- 
dices. Their basic graphs may be obtained from those of Group 1, 2 and 2', 3, 4, 7 and 10 
by making one simple change to them. For example, attaching the line G2 to quark line 
1 at a point below instead of above the line Gl on Group 2 and 2', or moving line G3 in 
Group 1 and 3 to above the gf/g- vertex on quark line 2, or flipping any of the basic graphs 
of Group 3, 4, 7 and 10 upside down. However, these extra groups may be included by 
giving a factor of two on each of their associated group because these are related to each 
other by a simple change of variables. Finally, we must mentioned that the graphs in 
each group must be subjected to further permutations of the three light quark lines to 
give further graphs. But these also can be taken care of by a numerical factor so each 
possible insertion point will generate only one graph and there is no need to divide the 
graphs any further or add any more groups. 

9 Colour Singlet And Octet Contribution In The Stan- 
dard Scheme 



9.1 Colour Singlet Contribution 

The procedure to obtain the colour singlet contribution in the standard scheme is quite 
similar to the one we used in Sec. |^. Remembering that in the standard scheme, the 
internal transverse momenta are taken to be negligible in comparison with the virtualities 
in all propagators in the perturbative hard part, so Th{x) is free from any fc^'s and the 
latter can be integrated over each wavefunction to give the corresponding distribution 
amplitude. So setting the fe^'s in ^j, Q and Qc in Eq. (pTl) to zero, is now independent 
of A;_L and kjj 

n(.,y) = 2V2 (4.)3 

9v3m5 xixsyiys (1 - xiy{l - yi) [xi{l - yi) + yi{l - xi)] 
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2xi) 



+ (x < — > y) 



2-J 



+ 



{-ly ^ xijxi - yi) 
-Vi) +1/1(1 -xi) 



(47) 



The renormalization scale fipt in the above equation has been set at the constant scale 
rric in because each gluon takes approximately ^ rric from the charmonium 

and so the virtuality is roughly m^. Now Fourier transforming the transverse momentum 
independent to transverse position-space yields several delta functions which force all 
6's and b"s to the origin. As a consequence, following from its definition and derivation, 
the Sudakov factor has to be set to unity. The decay form factor defined in Sec. | of xj 
into B — 13 becomes 



B{1) 



. v^|i?;(0)|aj 

ll 

'8v/?mc/^(47r)4 



[dx][dy] T'H^'\x,y) l|vl/^(x, 0)v[/^(y, 0) || (48) 



,B 



Bl 



where cxj = l/y^, 1 for J = 1,2 as before. This has to be combined with that of 
the colour octet contribution to be discussed below in accordance with our theoretical 
arguments given in Sec. || to give the true partial decay width. 



9.2 Colour Octet Contribution 

The total colour octet contribution to i3j has to be the sum over all contributing 
graphs from each group and over all possible helicity configurations of the light quarks 
and antiquarks in the outgoing baryon-antibaryon given in the appendices 



(49) 



Al,A2,A3=± 



The helicity dependent hard perturbative parts are 



{Ti(')(x,i/)}A„A,,A3 = ^(47ra,(m,))3V4rf^(2m,)' 



SgPgm{x,y){N^g^{x,y)}x^MM 

(50) 



' gm 

g^Groups 
M ember s 



where Pgm{x, y) is the product of propagators of the member m =U1, U2, U3, Ll, L2, L3, 
. . . etc. of the group g =1, 2, 2', . . ., 10. They can be obtained from the product of the 
basic propagators of each group as discussed in Sec. ||. The Sg is a symmetry factor for 
the group g to take care of similar potential groups that are related to g by simple change 
of variables. For the groups listed in the appendices, Sg ={2,2,2,4,2,1,1,2,1,1,2} for the 
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group g =1,2, 2', 3, . . ., 10, respectively. The coupling af"-^*, taken to be equal to vr, is that 
of attaching the constituent gluon to the basic graphs. It needs special treatment because 
of its different nature in comparison to the rest of the a^'s The {A^gm(x, ?/)}ai,a2,A3 



in the above equation contains the helicity dependence and denotes the numerator of the 
member m of the group g. These numerators are given in the appendices. Note that 
the calculation is really only gauge-invariant to order z"^ in the numerator (see ref. ||21|| ) 
because the cc are treated as on-shell, therefore all or higher terms in the numerator 
have to be dropped. The remaining helicity dependent product of the scalar functions 
^f23(x, 0) of the baryon wavefunctions are given below. 



||vi/^«(x,0)vi/^«(y,0)||+,+,_ =2 I ^fi,{x,0)^fi,{y,0) + ^g\{x,0)^g\{y,0) 

+{^fi2ix, 0) + ^g\ix, o)){^fUy, 0) + 0); 

'(51) 

||^^«(a:,0)^««(l/,0)||+,_,+ =2 | ^f,%ix,0)^fi,iy,0) + ^ii,ix,0)^i,\iy,0) 

+(^fl3(x, 0) + ^g\ix, o)){K%iy, 0) + ^iiiiy, 0)) } 

(52) 

||^^«(x,0)^^«(l/,0)||_,+,+ =2 { ^g%ix,0)^g%iy,0) + ^if,ix,0)^iUy,0) 

0) + ^iu^, o)){K%iy, 0) + ^iUy, o)) } 

(53) 



|v|>^(x, 0)^\y, 0)||+,+,_ = 6 I vj/A ^(x, 0)^i,{y, 0) + ^^x, 0)^i,{y, 0) 

+ {^,,{x, 0) - ^,,{x, 0))(^f32(2/, 0) - ^,,{y, 0) 



(54) 



The other two helicity arrangements of this quantity for the A follow the pattern above 
for the other octet baryons. For the decuplet baryons, they are 



|M/^^«(a;,0)vl/^^«(y,0)||+,+,_ = 3 | ^fi^{x,0)^fi^{y,0) + ^^i^{x,0)^^i^{y,0) 
I ^^1° {x, 0)^^^« {y, 0) II +,_,+ = 3 I ^f2^^ (x, 0)^f2^3° (y, 0) + i^, 0)^|\° iv, 0) 
|^^"'(x,0)^^"'(l/,0)||_,+,+ = 3 ( ^f/|(x,0)^f/|(y,0) + ^fi^^(x,0)^fi^2°(z/,0)^ . 



(55) 
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With these notations, the colour octet contribution looks rather simple. One only has to 
sum up all graphs from each and every group and then perform the integrations. 



10 The Widths Of xj Decay Into Baryons- Antibaryons 

With our method discussed in the previous sections, colour octet contribution which we 
argued to be necessary in additional to the singlet contribution to getting the correct 
P-wave Xj partial decay widths can be included. Before giving the numerical results, it 
must be mentioned that of all the baryons we considered, only the decay into proton- 
antiproton is measured so the majority of our results are in fact predictions. Moreover, 



the most recently reported measured values p9[ differ by as much as a factor of two from 



those in the Particle Data Tables |2^. This is also true in the case of the decay into 
pseudoscalar mesons as noted already in [|2T]]. Therefore we are content with results that 
lie somewhere in between these measurements. Until the situation improves, the colour 
octet decay constants fitted in [pO| , pT| cannot be more accurately determined hence also 



the current results. In any case, part of our goal is to show that explicit calculations do 
indeed support the theoretical arguments, that is the colour octet contribution in P-wave 
Xj decays cannot be neglected both in inclusive as well as in exclusive process. 

Our results are shown in Table |^ for octet and decuplet baryons. The numerical 
parameters used to obtain these results which were given throughout this paper have 
been collected together again in appendix |B| so that any interested readers do not have 
to search through the paper for their values. Only the kinematically plausible decuplet 
baryons with the lowest masses are considered. In getting this results, we used the colour 
octet decay constant f'^^ = 0.9 x lO^"' GeV^ fitted in [^] for the decay of X2 tttt within 
the standard hard scattering approach. We found that in order to obtain reasonable 
agreement of the Xi ~^ PP decay, it must have a smaller colour octet decay constant. 
Therefore we use fi^f = 0.225 X 10-3 GeV^. Note that this does not contradict the result 
in [^, since xi cannot decay into nn because of parity, and the decay constant fj^^ 



is therefore unconstrained in |2^. The smaller colour octet decay constant for the xi 
than those of the J = 0, 2 partners is not too surprising given that the odd-spin P-wave 
charmonium is a somewhat different heavy quark-antiquark system. Table ^ shows that 
the decay widths decrease roughly with increasing value of the baryon masses or with the 
reduction in the available phase space. 

In Table ^, the branching ratios of the decays into pp is shown against the experimental 
measurements. As mentioned above, the disagreement between PDG and BES is sizable. 
This is true not only in the ratios into pp but also in the case of a pair of pseudoscalars 
in the final state. Since our results are based on fit to the latters, until the experimental 
situation improves, we have to settle for our current results. 

From Sec. ^ and Table ^, it is clear that only by including colour octet can the decay 
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Octet 


r(i)+(8) (eV) 


Decuplet 


r(i)+(8) (eV) 


Baryons 


J=l 


J=2 


Baryons 


J=l 


J=2 


XJ^NN 


(56.27) 


154.19 


XJ AA 


33.49 


124.62 


xj 


28.42 


97.69 


XJ ^ 


18.46 


71.09 




21.49 


72.62 




9.42 


41.16 


Xj^ AA 


33.64 


69.19 









Table 3: The partial decay widths for xj decay into octet and decuplet baryon-antibaryon 
pairs. The width of xi — ^ in parenthesis is to indicate that this value is a fit unlike 
all partial widths of X2 which are predictions. Based on this fit, the rest of Xi widths are 
also predictions. 





Branching ratios (x 10 ^) 




PDG 


BES 


Xi PP 

X2 ^PP 


(6.39) 
7.71 


8.60 ± 1.2 
10.00 ± 1.0 


4.30 ± 2.3 ± 2.9 
5.90 ± 3.1 ± 3.3 



Table 4: Comparing our results with the measured widths from the PDG data and 
from the BES collaboration [^. This branching ratio of Xi is ^ fit- 



widths of the P-wave charmonia be brought in range of the experimental measurements. 
It would be better to show this in the modified hard scattering scheme which has the 
advantage of the dynamical setting of renormalization scales by the process itself. But 
due to the complexity of the calculations, we have to revert back to the standard scheme. 
Remembering this dynamical setting of scales is only possible if the Landau pole in as is 
suppressed by the Sudakov factor or better yet if it is not present at all. This suggests the 
use of one of the analytic models for which are free from the problem of the Landau 
pole |32|, 0]. The one suggested in |Q is particularly appealing because of its 
relative simplicity. By combining c^^naiytic ^j^]^ ^j^g standard hard scattering approach, a 
simpler scheme than the modified one can be constructed, but nevertheless still preserving 
the best features of that scheme. Under the new scheme, the amplitude would be given 

by 

M ~ fxAxA^) ® InMx) ® h^pNix) ® Th{x, ar'^'-ix)) . (56) 
We have also tried using this semi-modified scheme and again both colour singlet and 



octet must be included. The details will be given elsewhere |3^. A somewhat similar use 
of this analytic model in exclusive process to study pion form factors can be found in 
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Thus having satisfied ourselves with the genuine necessity and the correctness of the 
inclusion of the colour octet state in P-wave charmonium decays, we can now general- 
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ize the arguments to even higher wave quarkonia. Remembering that the P-wave QQ 
wavefunction has a 1/M suppression due to angular momentum, so one can deduce in a 
straight forward manner that for a D-wave, the wavefunction would be doubly suppressed 
by in relation to a S-wave. Then to calculate the exclusive decay of a D-wave 

quarkonium, one would need to include not only the colour singlet valence state and the 
next higher colour octet state now with the heavy fermions in a P-wave, but the next-next 
higher state must also be included for a consistent calculation. 

Finally as mentioned in the introduction, nucleon wavefunction constructed from QCD 
sum rules has large disagreement with the magnetic form factor measurements below 
50 GeV^. Therefore while lacking an alternative method to derive wavefunction in a 
fundamental way, we will have to satisfy with phenomenological constructions. We have 
shown together with fill] that the so-constructed wavefunction and its generalization to 
the flavour octet and decuplet baryons provide a set of reasonable model wavef unctions. 
They would no doubt provide a useful basis for future study of other exclusive processes 
involving baryons at moderate Q^. 
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Basic Propagators: 

1 1 1 1 1 

(2mc)2 X3(l-2/i)(2mc)2 xiyi{2mc)'^ (l-a;i)(l-2/i)(2mc)2 xsysi^rric)'^ 



Figure 7: Basic graph of Group 1 

A Appendix: Basic Graphs And Contributions For 
colour octet 

A.l Group 1 



LP. 


Num. for graphs of Group 1 with 
(Ai,A2,A3) = (+,+,-) 




J=2 


U3 








L3 


V2(x3 - z)(y3 - z){l + yi- z) 


U2 


-^xs{{x2 -z)(l + y,-z)- {-iy2zix, - y,)} 


M2 


-^xs{xs-z){l+yi-z) 


L2 


^{x3-z){y2-z){l + yi-z) 


Ul 


^X3{(xi -z)(l-yi- 2z) + (-iyz{xi - yi + 3(1 - z))} 


LI 


-^xsiil + yi - 2^)(yi - z) - (-1)^^(1 - 2yi + z)} 


G3 


-V2{x3 - y3){x3 - z){l + yi- z) 


G2 


^X3{ 2(1 - z^) - yi(3 + yi - &z) + a;i(3 - yi - 2z) 
-(-l)^2(xi-yi)(l-yi-2^)} 


Gl 


^a;3{2(a;i-yi)(l+yi-2z) 

+(-l)-^(2^(2 - ^) + a;i(l - yi + 2z) + yi(3 - yi - 6^))} 


4G 





l(iv'2 
, .'3 



Table 5: (a) Numerators of the graphs of group 1 with hehcity (+,+,-). 
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LP. 


Num. for graphs of Group 1 with 
(Ai,A2,A3) = (+,-,+) 


J=l 


J=2 


u o 








T 

J_iO 


(+ + -) + (- + +) 


U2 


-^xs{x2-z){l+y,-z) 


M2 


-z){l + y,-z) + {-iy2{,; y,){l -y,- z)} 


L2 


^{x3-z){y2-z){l + yi-z) 


Ul 


^xs{{x, -z){l-y,- 2z) + {-iyz{x, -y, + 3(1 - z))} 


LI 


-^xsiil + yi - 2z){y, - z) - (-l)'^^(l - 2y, + z)} 


G3 


(+,+,-) + (-,+,+) 


G2 


^X3{ 2(1 - z^) -y,{3 + yi- 6z) + xi(3 - t/i - 2z) 
-{-iy2{x^-y,){l-y,-2z)} 


Gl 


^xs{2{x,-y,){l+y,-2z) 

+(-l)-^(2z(2 -z)+ xi{l - yi + 2z) + yi{3 - yi - 6z))} 


4G 





8 ^3 



Table 5: (b) Numerators of the graphs of group 1 with hehcity (+,-,+). 



LP. 


Num. for graphs of Group 1 with 

(Ai,A2,A3) = (-,+,+) 


X-Prop. 


J=l 


J=2 


U3 








n. g. 


L3 


-{-iy2V2z{x,-y,){l-y,-z) 


1 

-2(j/3-2)(2m,c)2 


U2 








—z(x2~z)(2m,c)'^ 


M2 








X3(l-S/l)(2me)2 


L2 


(_l)J2^;,(a;,_^^)(l_^^_^) 


1 


Ul 








—z(xi — z)(2mc)'^ 


LI 








-z(jn-z)(2rnc)^ 


G3 


-(-l)-^v/2(,r,-y0 


{.r, + z){l-y,-z) 


1 


G2 








1 

(l-xi)(l-2/i)(2mc)^ 


Gl 








xiyi (2mc)^ 


4G 








1 

(2mc)2 



Table 5: (c) Numerators and the additional propagators of the graphs of group 1 with 
hehcity (-,+,+). 
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Basic Propagators: 

11 1 11 

(2mc)2 (l-2/i)(2mc)2 (l-j;i)(l-j/i)(2mc)2 X3(l-yi){2mc)'^ xsysiirric)'^ 




SiSL 





Figure 8: Basic graph of Group 2 



A. 2 Group 2 and Group 2' 
A.2.1 Group 2 



LP. 


Num. for graphs of Group 2 with (Ai, A2, A3) 


(+,+,-) 


(+,-,+) 


J=l 


J=2 


J=l 


J=2 


U3 














L3 








( ly^lfzil y, zf 


U2 










M2 








( iy'\fx,{l y, zf 


L2 











Ul 


(-iy^-^x,z{l-y,-z) 


Ml 














LI 














G3 














G2 


'-^xs{{l-zmi-y,)-z) 

-{-iy{l-yi-2z){l-y,-z)} 



Table 6: (a) Numerators of the graphs of group 2 with hehcity (+,+,-) and (+,-,+). 
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LP. 


Num. for graphs of Group 2 with 

(Ai,A2,A3) = (-,+,+) 


X-Prop. 


J=l 


J=2 


U3 








n. g. 


L3 




1 

— z{y;i, — z){2mcy^ 


U2 








— z{x2 — z){2inc)'^ 


M2 








x-3(l-i/i)(2mc)2 


L2 


( l)^32,f ;.(! y, zf 


1 

-z(y2-z){2mc)'^ 


Ul 








1 

—z(xi—z)(2mc)'^ 


Ml 








n. g. 


LI 








n. g. 


G3 








n. g. 


G2 








1 

(l-j;i)(l-j/i)(2mc)2 



Table 6: (b) Numerators and the additional propagators of the graphs of group 2 with 
hehcity (-,+,+). 
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Basic Propagators: 












11 1 11 

(2mc)2 (l-2/i)(2mc)2 (l-j;i)(l-j/i)(2mc)2 {l-xi)y3{2mc)'^ xsysiirric)'^ 



Figure 9: Basic graph of Group 2' 



A.2.2 Group 2' 



LP. 


Num. for graphs of Group 2' with 
(Ai,A2,A3) = (+,+,-) 




J=l J=2 


U3 




L3 




U2 


-l^(l-z)(x2-z){ys-z) 


M2 


-^^y,{l-z){y,^z) 


L2 


-'-^ysiil - z){y, -z) + (-l)^z(l - 1/1 - z)} 


Ul 


-'-^y,{l-z){x,-z) 


Ml 


^y,{{l-z)-{-iy{l-y^-z)} 


LI 




G3 


2y/2{l -z){x^-y^){yz-z) 


G2 


-^yz{{l-z){x,-y,-?>z) 

-{-iy{l-x^-2z){l-y^-z)] 



Table 6: (a') Numerators of the graphs of group 2' with hehcity (+,+,-). 
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LP. 


Num. for graphs of Group 2' with 




(Ai,A2,A3) = (+,-,+) 




J=l J=2 


U3 


(+,+,-) + (-,+,+) 


L3 


-^J^H - z){y, - zf 


U2 


-^-^{l-z){x,-z){y,-z) 


M2 


-^^y,{{l-z){y,-z) 




-{-iy{l-x^-z){l-y^-z)} 


L2 


7™ 

-^J^y,{l-z){y,-z) 


Ul 


-^-^y,{l - Z){X, - Z) 


Ml 


^-^y,{{l-z)-{-iy{l-y,-z)} 


LI 


'-^y,{y,-z + {-iyz} 


G3 


(+,+,-) + (+,-,+) 


G2 


-^yz{{l-z){xi-yi-'iz) 




-{-iy{l-x^-2z){l-y^-z)] 



Table 6: (b') Numerators of the graphs of group 2' with helicity (+,-,+). 



LP. 


Num. for graphs of Group 2' with 
(Ai,A2,A3) = (-,+,+) 


X-Prop. 


J=l J=2 


U3 


{-iy2^z{l-x,-z){l-y,-z) 


1 

—2(0:3— z)(2mc)^ 


L3 








-z(2/3-z)(2mc)^ 


U2 


{-iyi^z{l-x,-z){l-y,-z) 


1 

—z(x2—z)(2mc)'^ 


M2 








(l-a:i)j/3(2mc)2 


L2 










Ul 








—z{xi—z){2mc)'^ 


Ml 








(l-tjl)(2in, 


LI 








1 

-z{y\-z)(2mc)'^ 


G3 


(-l)-^y2(l - - z){l - 1/1 - z){y^ + z) 


1 

xzyz(2mc)'^ 


G2 








(l-a:i)(l-j/i)(2mc)2 



Table 6: (c') Numerators and the additional propagators of the graphs of group 2' with 
helicity (-,+,+). 
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Basic Propagators: 



1 1 1 

{(2i-a;i){zi-2/i)-l/4}(2mc)2 X3(l-yi){2mc)'^ xiyi(2mc)^ 



{l-xi)(l-yi){2mc)^ X3y3{2mcY 



Figure 10: Basic graph of Group 3 
A.3 Group 3 



LP. 


Num. for graphs of Group 3 with 
(Ai,A2,A3) = (+,+,-) 


J=l J=2 


U3 




L3 




U2 


-^x^{{x2 - z){2y^ + z)- {-iy2z{x^ - y^)} 


M2 


-^x^{x^-z){2y^ + z) 


L2 


-^{x;~z){y2-z){2y, + z) 


Ul 


^^x^{{x,-z){2y,-z) + {-iy2z] 


LI 


-^x^{{2y, - z){y, - z) - {-iy2z{l - 2y, + z)} 


03 




G2 


-^^3{ {xi +yi- 2z)i2y, + z)- 2(2x1 + z) 

+{-iy{xi-yi){l-yi-2z)} 


Gl 


^X3{{x,-y,)i2y,-z) 

-i-iyHxi + z){2y, -z)- {x, +1/1 + 2;^))} 


Q 


-^x^i \{2yi + z){2yi - z) - [x^ ~Vyi^z) 

+ (-1)^(1(2x1 + z){2y^ -z)- (xi + yi + z))} 



Table 7: (a) Numerators of the graphs of group 3 with hehcity (+,+,-). 
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LP. 


Num. for graphs of Group 3 with 
(Ai,A2,A3) = (+,-,+) 


J=l J=2 


U3 




L3 


^{(X3 - z){y^ - z){2y, + - (-l)^2z(xi - y,){l -y,- z)} 


U2 


-^xs{x2-z){2y, + z) 


M2 


-^x^Uxs - z)(2yi + z) + (-iy2(xi - yi)(l - yi - z)} 


L2 


-^ix3-z)iy2-z){2y, + z) 


III 


''^'^T-.liT^ z)(2v, z) + ( iy2z\ 


LI 


-^x^{{2y, - z){y, - z) - (-1)'^2^(1 - 2y, + z)} 


G3 


(+, +, -) + +• +) 


G2 


-^X3{ (,xi + - 2;2)(2|/i + z)- 2(2x1 + ;2) 
+(-l)^(xi-|/i)(l-yi-2^)} 


Gl 


^x;{{x,~yi){2yi-z) 

-(-l)''((.ri + .:)(2yi---)-(.ri + yi + 2,:))} 


Q 


-||x3{ |(2|/i + z){2y^ -z)- (xi + 1/1 + ^) 

+(-l)-^(i(2xi + z){2y^ -z)- (xi + yi + z))} 



Table 7: (b) Numerators of the graphs of group 3 with hehcity (+,-,+). 



LP. 


Num. for graphs of Group 3 with 

(Ai,A2,A3) = (-,+,+) 


X-Prop. 


J=l 


J=2 


U3 








1 

—z{x3—z)(2mc)'^ 


L3 


-{-lYl^z{x,-y^){l-y,-z) 


1 

-z(y3-z)(2mcY 


U2 








—z{x2 — z){2mc)^ 


M2 








X3(l-J/l)(2me)2 


L2 


-{-lY^-^z{x,-y,){l-y,-z) 


1 

-2(j/2-2;)(2mc)^ 


Ul 








— z{xi-z){2mc)'^ 


LI 








-z(y\-z)(2mc)'^ 


G3 


-(-1)^^(X3 + z){x, - y,){l -yi-z) 


1 

X3ys,(2mcY 


G2 








(l-a:i)(l-j/i)(2mc)2 


Gl 








xiyi(2mcY 


Q 








[(xi-2i)toi-2i)-l/4l(2mc)2 



Table 7: (c) Numerators and the additional propagators of the graphs of group 3 with 
hehcity (-,+,+). 
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Basic Propagators: 



1 1 1 

{(22-a;2){22-2/2)-l/4}{2mc)2 (l-a;2)(l-2/2)(2mc)2+p2 xiyi{2mc)'^ 
1 1 



X 



X2j/2(2mc)^ X3y3(2mc)'^ 



Figure 11: Basic graph of Group 4 



A.4 Group 4 



LP. 


Num. for graphs of Group 4 with (Ai, A2, A3) 


(+,+,-) 


(+,-,+) 


J=l J=2 


J=l J=2 


U3 


-^{X3- Z){ (.Ti - ,T3 + Z) + {yi - ?/3 + Z) 

-2{xi -X3 + z){z2 - 1/2) 
-2(yi - + z){z2 - X2)} 


^-iy^z{x2-y2){2x^ + x^-z) 


L3 


^{ys - z){ {xi -X3 + z) + {yi -y3 + z) 
-2{xi -X3 + z){z2 - y2) 
-2(yi -y3 + z){z2 -X2)} 


(-l)-^f z(x2-y2)(2yi + y3-^) 


U2 














L2 














Ul 


(+,-,+) + (-,+,+) 


{^_lY^z[x2-y2){xi + 2x^-z) 


LI 


(+,-,+) + (-,+,+) 


^_iy:^z{x2-y2){yi + '2yz-z) 


G3 


- y3){ {xi - X3 + z) + (yi -yz + z) 
-2{xi - x^ + z){z2 - y2) 
-2(yi -y3 + ^)(^2 -2:2)} 


(-I)^^(x2-y2) 

X {{2xi + x^ - z){y^ + z) 
{2yr + y:,-z){x:, + z)} 


G2 








Gl 


(+,-,+) + (-,+,+) 


{-iy^-§{x2-y2) 
X {(2a;3 + a;i -^)(|/i + 2;) 

{^yz + yi - z){xi + z)} 


GR 














4G 


-^{x2-y2) 


-^(^2-^2) 


-(-l)-^¥(x2-y2) 


Q 












Table 8: (a) Numerators of the graphs of group 4 with hehcity (+,+,-) and (+,-,+). 



37 



LP. 


Num. for graphs of Group 4 with 
(Ai,A2,A3) = (-,+,+) 


X-Prop. 


J=l J=2 


U3 


(+,+,-) + (+,-,+) 


1 


L3 


(+,+,-) + (+,-,+) 


-z(y3-z)(2mc)'^ 


U2 








n. g. 


L2 








n. g. 


Ul 


-^(xi - z){ (xg -xi + z) + (ya -yi + z) 
-2{x3 -xi + z){z2 - 1/2) 

-2(1/3 - 1/1 + ^)(^2 - 2:2)} 


1 

— z)(2mc)^ 


LI 


- ^){ (a^a - xi + z) + (t/3 - t/i + 2;) 

-2(,T3 - Xi + 2;) (2:2 - 7/2) 

-2(y3-yi + -)(.:2-.r2)} 


1 


G3 


(+,+,-) + (+,-,+) 


1 

0:32/3 (2mc)2 


G2 





n. g. 


Gl 


-^{xi - yi){ {x3 -xi + z) + (1/3 -yi + z) 
-2(2:3 - a;i + z){z2 - 1/2) 
-2(l/3 - 1/1 + z){z2 - X2)} 


1 


GR 








n. g. 


4G 


-73 (.'•2 - <J2) 




1 


Q 








n. g. 



Table 8: (b) Numerators and the additional propagators of the graphs of group 4 with 
hehcity (-,+,+). 
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Basic Propagators: 



1 



1 



{(2i-a;i){zi-2/i)-l/4}(2mc)2 {(z2-a;3)(z2-J/3)-l/4}(2mc)2 xiyi{2mc)^ 
1 1 



X 



X2j/2(2mc)^ X3y3(2mc)'^ 



Figure 12: Basic graph of Group 5 



A. 5 Group 5 



LP. 


Num. for graphs of Group 5 with 
(Ai,A2,A3) = (+,+,-) 


J=l 


J=2 


U3 








L3 








U2 


-{-iy^z{2xi + z){2x3 + z) 


L2 


-{-iy^zi2y, + z)i2ys + z) 


Ul 




LI 


-(-1)-^^^{1 - 2^2 - 4(1/2 + ys - z,){xs - z,)} 


G3 








G2 








Gl 


+2X2(21/3 - 

-2y2(2x3 - 


I - 4(X3 - Zi)(2/3 - 2:1)) 

^z){yi + z) 
^z){xi + z)} 


UQ 


(^3 - Z/3) - M^2 - Z/2) + 4(X1 - yi)(l - Z,f 

-4(xi(l - X3) - X22;i)(l - yi - zi) 
+4(1/1(1 - ys) - y2Zi){l -xi- zi)} 


LQ 









Table 9: (a) Numerators of the graphs of group 5 with hehcity (+,+,-). 
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LP. 


Num. for graphs of Group 5 with 
(Ai,A2,A3) = (+,-,+) 




J=l 


J=2 


U3 


-^{ 2(a;3 -z)[l-X2-y2- 2(1 -x^- zi){yi - z^) 
- 2(1 -ys- zi){xi - zi)] 
+{-iyz[l - 2x2 - 4(1 - X3 - 2i)(yi - zi)]} 


L3 


^{ 2(1/3 -z)[l-X2-y2- 2(1 - X3 - z,){yi - Zi) 
- 2(1 - 1/3 - Zi){xi - Zi)] 

+ (-l)-^^[l - 2y2 - 4(1 - 1/3 - z,){x, - z,)]} 


U2 


-f{x2-z){l- {x2 + y2-2z) 
+2{x3 - zi){yi - zi) +2{ys - zi){xi - zi)} 


L2 


-f{y^-z){l- {x2 + y2-2z) 
+2{x3 - zi){yi - zi) +2(i/3 - zi){xi - zi)} 


Ul 


^{2{xi - z)[l - X2 - 1/2 - 2(1 - xi - z^){ys - z^) 

- 2(1 - 1/1 - Zi){x3 - Zi)] 
+ (-l)'^^[l - 2,X2 - 4(1 -X,- z^){y3 - zi)]} 


LI 


-^{ 2(1/1 -z)[l- X2 - 1/2 - 2(1 - xi - z^){y3 - z,) 
- 2(1 - yi - Zi){x3 - Zi)] 
+ (-l)-^2[l - 2y2 - 4(1 - yi - 2i)(x3 - ^i)]} 


G3 


4(x3 - 1/3) [1 - 1/2 - X2{2yi + z)- 2{xi - zi){yi - ya + z)] 
-{-iy[2z{x2{2yi + z)- 1/2(2x1 + z)) 
+0:3(1 - 2y2 - 4:{xi - zi){l - 1/3 - zi)) 
-ysil - 2x2 - 4(yi - zi){l - xs - zi))]} 


G2 








Gl 


^{4(a;i - yi)[l - 1/2 - 2:2(21/3 + z) + 2(2:3 - Zi)(yi - y^ - z)] 

^{-iy[2z{x2{2y3 + z)- 1/2(2x3 + z)) 
+xi(l - 2i/2 - 4(a:3 - 2;i)(l - yi - Zi)) 
-1/1(1 - 2x2 - 4(1/3 - Zi){l -Xi- Zi))]} 


UQ 


-^{ z{2 + z){xi - yi) + {z{2 -z)+ 2(xi + yi))(x2 - ^2) 
+4(zi - y3)xl - 4(zi - X3)yf } 


LQ 


^{ z(2 + z)(x3 - ys) + (z(2 -z) + 2(xs + y3)){x2 - ^2) 
+4(^1 - yi)xl - 4(zi - xi)y|} 



Table 9: (b) Numerators of the graphs of group 5 with hehcity (+,-,+). 
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LP. 


Num. for graphs of Group 5 with 
(Ai,A2,A3) = (-,+,+) 


X-Prop. 


J=l J=2 


U3 


-{-IY^Z{1 - 2X2 - 4(xi +X2- - Zl)} 


1 

-z{xz-z)(2mcY 


L3 


{-iy^z{l - 2|/2 - 4(1/1 + 1/2 - zr){x^ - z^)} 


1 

-2(j/3-z)(2mc)2 


U2 


-{-iy^z{2x^ + z){2x^ + z) 


1 

—z(x2—z)(2ma)'^ 


L2 


-{-iy^z{2yr + z){2y^ + z) 


1 

-z{y2-z)(2mcY 


Ul 








—z(xi—z){2mc)'^ 


LI 








-z(yi-z)(2mcY 


G3 


(-1) (^3 - y,){l - 4(a:i - z,){y, - z,)) 
+2x2(2yi + z)(y3 + z) 
-2y2{2xi + z){xz + z)} 


X3yi{2mcY 


G2 








n. g. 


Gl 








1 

a;ii/i(2mc)2 


UQ 








Ka:i-2i)(yi-^i)-l/41(2mc)2 


LQ 


(^1 - y^) - 2^(^2 - y2) + 4(X3 - y^yi - z^f 

-4(,T3(1 - Xi) - X2Zi){l - 7/3 - Zi) 

+ i(y.s(i-yi)-y2--i)(i-.rs - --i)} 


1 

[{a:3-zi)(2/3-zi)-l/4](2mc)2 



Table 9: (c) Numerators and the additional propagators of the graphs of group 5 with 
hehcity (-,+,+). 
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Basic Propagators: 

1 1 1 1 

(2mc)* xiyi{2mc)'^ 122/2 (Smc)^ a;32/3(2mc)^ 



Figure 13: Basic graph of Group 6 
A. 6 Group 6 



LP. 


Num. for graphs of Group 6 with (Ai, A2, A3) 


(+,+,-) 


(+,-,+) 


J=l J=2 


J=l J=2 


U3 


2^2(0:3 - z) 


-V2{xs - z + {-iy2z} 


L3 


2v^(2/3 - z) 


-V2{y3-z + i~iy2z} 


U2 


-V2{x2-z+{-iy2z} 


2V2{x2 - z) 


L2 


-^2{!j,-z + {-iy2z} 


2/2(^2 -.:) 


Ul 


-V2{xi-z+{-iy2z} 


LI 


-V2{yi-z + {-iy2z} 


G3 














G2 














Gl 















Table 10: (a) Numerators of the graphs of group 6 with hehcity (+,+,-) and (+,-,+). 
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LP. 


Num. for graphs of Group 6 with 
(Ai,A2,A3) = (-,+,+) 


X-Prop. 


J=l J=2 


U3 


-V2{x3- z+{-iy2z} 


1 

—z(xz—z)(2mcY 


L3 


-V2{y3-z + i-iy2z} 


1 

-z(y3-z)(2mc)^ 


U2 


-V2{x2 - z + {-iy2z} 


1 

—z{x'2,—z){2incY 


L2 


-V2{y2 - z + {-iy2z} 


1 

-^(2/2-z)(2mc)^ 


Ul 


2\f2{xx - z) 


1 

— z(a;i— z)(2mc)^ 


LI 


2^2(1/1 - z) 


1 

-z{y\—z){2mc)'^ 


G3 








n. g. 


G2 








n. g. 


Gl 








n. g. 



Table 10: (b) Numerators and the additional propagators of the graphs of group 6 with 
helicity (-,+,+). 
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Basic Propagators: 

11 1 11 

(2mc)2 (l-2/i)(2mc)2 (l-X2)2/3(2mc)2 X2y2(2mc)2 X32/3(2mc)2 









Figure 14: Basic graph of Group 7 



A. 7 Group 7 



LP. 


Num. for graphs of Group 7 with (Ai, A2, A3) 




(+,+,-) 


(+,- 


-,+) 




J=l J=2 


J=l 


J=2 


U3 




-i-iy'-^z 


(1 - X2 - z) 






x{l-yi-z) 


L3 


^-^{1 - z){y, - zf 








U2 


-^^{l-z){x,-z)ys 








L2 


'^y3{{l-z){y2-z) 

+(-iyz{l-y,-z)} 








Ul 


^-^{l-z){x,-z){y,-z) 


( l)'^26V2^(i ^) 






X 




UMl 


-^-^y,(l-z)iys-z) 








LMl 


-'-^ys{il-z)-i-iy{l-y,-z)} 








LI 


-'-^ys{{yi-z) + i-iyz} 








G3 


-^{l-z){xs-ys){ys-z) 


(-l)'^f (1 


- X2- z) 






x(l -1/1 -2;) 
x(y3 + ^) 


G2 


-ify,{2{l-z){x,-y,) 

-{-iy{x2 + z){l-yi-z)} 









Table 11: (a) Numerators of the graphs of group 7 with hehcity (+,+,-) and (+,-,+). 
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LP. 


Num. for graphs of Group 7 with 
(Ai,A2,A3) = (-,+,+) 


X-Prop. 


J=l J=2 


U3 


-[-lYl^z{l-X2-z){l-y,-z) 


1 

-z(x:i-z)(2mc)'^ 


L3 








1 

-z(y3-z)(2mc)'-^ 


U2 








—z(x2—z){2mc)^ 


L2 


(-1)^^^1/3(1 -yi-^) 


1 

—z{y2—z){2mc)^ 


Ul 








— z(x\ —z)(2mc)^ 


UMl 




1 

(l~X2)y3.(2mcY 


LMl 








(l-yi)(2mc)2 


LI 








-z(?/i-z)(2mc)^ 


G3 


(-l)''f (l-.r2-c)(l-//i-0(y. + .) 


1 


G2 


(-l)^^Z/3(l-Z/l-^)(^2 + ^) 


1 

X2y2(2mc)'^ 



Table 11: (b) Numerators and the additional propagators of the graphs of group 7 with 
hehcity (-,+,+). 



45 







Basic Propagators: 

11 111 

(2mc)2 a;2(l-2/3)(2mc)^ (l-3;2)2/3(2mc)2 a;22/2(2mc)^ a;32/3{2mc)2 









Figure 15: Basic graph of Group 8 



A. 8 Group 8 



LP. 


Num. for graphs of Group 8 with (Ai, A2, A3) 


(+,+,-) 


(+,-,+) 


J=l 


J=2 


J=l J=2 


U3 








( ly^'^i^X^zil X2 Z) 


L3 














U2 














L2 










Ul 








{-iy^,;,{:i-,;-y 


UMl 














LMl 














LI 










G3 








(-iy^-^x,{l-X2-z){y, + z) 


G2 


{-iy*-^y,(l-y,-z){x2+z) 









Table 12: (a) Numerators of the graphs of group 8 with hehcity (+,+,-) and (+,-,+). 
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LP. 


Num. for graphs of Group 8 with 
(Ai,A2,A3) = (-,+,+) 


X-Prop. 

IT 


J=l J=2 


U3 


^X2{{xs - z){ys -z) + {-iyz{l -X2- z)} 


1 


L3 


-'-^x^iys - zf 


1 

-2(j/3-z)(2mc)^ 


U2 


-'^y,{x, - zf 


1 

-z{x2-z){2mc)'^ 


L2 


^yA{x2 - z){y2 -z) + {-iyz{l - 1/3 - z)} 


1 

-z(:y2-2)(2mc)2 


Ul 


-^X2{xi - z)iy3 - z) 


1 

—z(xi—z)(2mc)''^ 


UMl 


-^X2|/3{(l/3 -Z)- (-l)-^(l - X2 - Z)} 


1 

(l-a;2)j/3(2mc)2 


LMl 


-'-^X2ys{{x2 -Z)- {-iy{l - 1/3 - Z)} 


1 


LI 


-^y3{x2 - z){y, - z) 


1 

-z{yi-z){2mc)'^ 


G3 


^X2{2{x, -y,){y,-z) 

+ {-iy{l-X2-z){ys + z)} 


1 


G2 


-^y,{2{x2-y2){x2-z) 

-(-i)-^(i-y.---)(.r2 + .--)} 


1 

X'zyii'ii'ic)'^ 



Table 12: (b) Numerators and the additional propagators of the graphs of group 8 with 
hehcity (-,+,+). 
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Basic Propagators: 



(2mc)2 (l-a;2)(l-j/2){2mc)2+p2 xiyi(2mc)'^ X22/2(2mc)2 



X3?/3(2mc)^ 



Figure 16: Basic graph of Group 9 

A.9 Group 9 



LP. 


Num. 


for graphs of Group 9 with (Ai, A2, A3) 






(+,+ 


-) 


(+,- 


,+) 




J=l 


J=2 


J=l 


J=2 


U3 


"75(^3 -^){(1+ 3^2 - 


- z){yi -ys + z) 


{-iyV2z{x2 - 2/2) (2x1 + X3-Z) 




+ (1 + 2/2 - z){xi - X3 + z)} 






L3 


-73(^3 -^){ (1+2:2 - 


-z){yi -ys + z) 


{-iyV2z{y2-X2){2y, + ys-z) 




+(1 + 2/2 - z){xi -X3 + z)} 






U2 














L2 














Ul 


-V2{ i(a;i - z)[{l + X2- z)(jj3 -yi + z) 
+(1 + 2/2 - 2;)(a;3 - xi + 2;)] 
-(-l)"'2;(a;2 - 2/2)(2a;3 + xi - ^)} 


(-1)'^V2^(X2 - 2/2) (2X3 + xi-z) 


LI 


-V2{ 1(2/1 - z)[(l + X2 - ^)(2/3 -yi + z) 

+ (1 + 2/2 - 2;) (2:3 - + z)] 
-{-iyz{y2-X2){2y3 + yi-z)} 


{-iyV2z{y2-X2){2y3 + yi-z) 


G3 














G2 














Gl 














GR 














4G1 














4G2 















Table 13: (a) Numerators of the graphs of group 9 with hehcity (+,+,-) and (+,-,+). 
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T P 


Num. for graphs of Group 9 with 

(Ai,A2,A3) = (-,+,+) 


A.-r rop. 


J=l J=2 




-V2{ ^{xs - z)[{l + X2- z){yi -ys + z) 

+ (1 + 1/2 - Z){xi -X3 + Z)] 

-{-iyz{x2 - 2/2)(2a;i + xs - z)} 


1 

— z(x3— z)(2mc)^ 


L3 


-V2{^iy3-z)[{l + X2-z){yi-ys + z) 

+ (1 +y2- Z){xi - ,T3 + Z)] 

-{-iyz{y2 - X2){2yi + 1/3 - z)} 


1 

-z{y3-z)(2incY 


U2 








n. g. 


L2 








n. g. 


Ul 


-^{xi - z){ {I + X2- z){yz -yi + z) 

+(1 + y2- z){x3 -xi + z)} 


1 

—z{x\- z)(2mc)'^ 


LI 


"73(^1 ~ (1 + ^2 - z){y3 -yi + z) 

+(1 + y2- Z){X3 -xi + z)} 


1 

~z{yi-z){2mc)'^ 


G3 








n. g. 


G2 








n. g. 


Gl 








n. g. 


GR 








n. g. 


4G1 








none 


4G2 








none 



Table 13: (b) Numerators and the additional propagators of the graphs of group 9 with 
helicity (-,+,+). 
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Basic Propagators: 



(2mc)2 (l-2/i)(2mc)2 {l-xi){l-yi){2mc)'^+p^ X22/2(2mc)2 



X3?/3(2mc)^ 



Figure 17: Basic graph of Group 10 

A.IO Group 10 



T T ) 

l.r. 


Num. for graphs of Group 10 with (Ai, A2, A3) 


(+,+,-) 


(+,-,+) 


J=l J=2 


J=l J=2 


U3 


\/2(l - z){x3 - z){y3 -y2- z) 


y^{(l-z)ix3-z){y3-y2-z) 

+ {-iyz{y2 + 2/3 - ^)(2a;2 + X3 - z)} 


L3 








U2 


V2{{l-z){x2-z)(y2-y3-z) 
+{-iyz{y2 + y3- z)(2x3 + X2- z)} 


V2{1 - z){x2 - z){y2 -y3- z) 


L2 














Ul 














Ml 














LI 














G3 


V2{1 - z){x3 - y3){y3 -y2- z) 


V2{(1 - Z){X3 - 2/3)(2/3 - 2/2 - ^) 

+ (-1)^(2/2 + 2/3-^) 

X {x3y2 + X2y3 + 2:32/3 

+ z(x2+2/2-^))} 


G2 


\/2{{l - z){x2 - 2/2) (2/2 -ys- z) 
+ {-iy{y2 + y3-z) 

X {X3y2 + X2y3 + X2y2 

+ z{x3 + y3 - z))} 


V2(l - 2;)(a;2 - 2/2) (2/2 - 2/3 - ^) 


GR 











4G 


V2{2{l-z)-{-iy{y2 + y3-z)} 


V^{2(l-2;)~(-l)-^(2/2 + 2/3-^)} 



Table 14: (a) Numerators of the graphs of group 10 with hehcity (+,+,-) and (+,-,+). 



50 



LP. 


Num. for graphs of Group 10 with 
(Ai,A2,A3) = (-,+,+) 


X-Prop. 


J=l 




U3 


{-iyV2z{y2 + 2/3 - z){2x2 + xs-z) 


1 

-z(xz-z)(2mcY 


L3 








n. g. 


U2 


{-iyV2z{y2 + 2/3 - z){2x3 + X2 - z) 


1 

-z(2/3-2;)(2mc)2 


L2 








n. g. 


TT1 

U 1 








n. g. 


Ml 

iVl J. 








n. g. 


LI 








11. g. 


G3 


{-iyV2{y2 + y3-z) 
X {x3y2 + X2y3 + X3y3 + z{x2 + y2 - z)) 


1 

X3yz{2mc)'^ 


G2 


{-iyV2{y2 + y3-z) 
X {x3y2 + X2y3 + X2y2 + z{x3 + y3 - z)) 


1 

X2y2{2mcY 


GR 





n. g. 


4G 


(-iy2V2{y2 + y3-z) 


1 

(2mc)2 



Table 14: (b) Numerators and the additional propagators of the graphs of group 10 with 
hehcity (-,+,+). 
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B Numerical Parameters 

Here we gather together our input as well as other parameters used in the calculations in 
table form. 

B.l Basic Parameters 



Symbol 


Value 


Aqcd 


0.22 GeV 




1.50 GeV 




1.00 GeV 


Hr (SHSA) 


mc 


^iR (MHSA) 


square root of the largest 




virtuality, see section ^ 


^/ 


4 



Table 15: Basic input parameters. 



B.2 Charmonium Parameters 



Symbol 


Value 


Origin 


\R'pm 


0.220 GeV^/^ 


see ref. S P 




0.225 X 10^^ GeV^ 


obtained from fit here 


7^^^ 


0.900 X 10-3 GeV^ 


from ref. 


m 


z = 


0.150 


from ref. 




Zi^ Z2 = {I - z)/2 


0.425 


from ref. 


10,11 



Table 16: Charmonium parameters. 
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B.3 Baryon Wavef unction Parameters 



Baryon 


5i 


B2 


^3 


B, 


^5 




0.750 


0.250 


0.000 


0.000 


0.000 


S 


0.216 


0.394 


-0.293 


-0.914 


0.241 




1.106 


0.050 


-0.282 


1.717 


-0.498 


A 


-0.721 


0.389 


-0.150 


-0.574 


0.093 


A 


0.000 


0.000 


0.000 


0.000 


0.000 


S* 


-0.547 


0.182 


-0.216 


-1.081 


0.062 




0.540 


-0.180 


-0.382 


1.742 


-0.413 



Table 17: Baryon wavefunction parameters derived in ref. with the constituent 

strange quark mass = 350 MeV, the octet baryon decay constant fs^ = 6.64 x 10~^ 
GeV^ and the transverse size parameter a^g = 0.75 GeV~^ at the reference scale /iq. The 
same for the decuplet baryons are fs^o = 0.0143 GeV^ and a^^o = 0.80 GeV~^. Note that 
in ref. |TT]] of all the decuplet baryons, only the parameters of A were given. 



From these parameters, the mean-squared internal transverse momentum of the baryons 
= (fc±^) can be worked out. The average value is p(8) = 415.0 MeV for the octet and 
P(io) = 389.0 MeV for the decuplet baryons. Their use as infrared cutoff was discussed in 
section |. 
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Figure Captions 

0) Basic graphs that could contribute to xj colour singlet decay. But actually, only 
graphs of type (a) can contribute. 

Graphs of type Fig. [I| (a) can be divided further into four groups. 

In addition to the graphs of type Fig. |l], these form the bases of further contributions 
in the colour octet decay channel. 

Our labelling scheme as applied to (a) Group 2 and (b) Group 4. 

Examples of complete colour octet graphs, (a) graph Ul and (b) graph L2 of Group 
5, and (c) graph L3 of Group 1. 

Examples of complete colour octet graphs with 4-gluon vertex, (a) graph 4G of 
Group 4 and (b) graph 4G of Group 10. 
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Table Captions 

The expansion coefficients of the distribution amphtudes (1)^23 of the octet baryons 
considered in the xj decay. The parameters associated with this set of coefficents 
are fsuM = 0.0143 GeV^ and ae.o = 0-80 GeV'^ 

Clearly, the colour singlet contributions are insufficient in explaining the experimen- 
tal data of Xj decay into pp. 

The partial decay widths for xj decay into octet and decuplet baryon-antibaryon 
pairs. The width of xi — ^ in parenthesis is to indicate that this value is a fit 
unlike all partial widths of X2 which are predictions. Based on this fit, the rest of 
Xi widths are also predictions. 



Comparing our results with the measured widths from the PDG |£3l data and from 
the BES collaboration E^. This branching ratio of Xi is a fit. 



Numerators of the graphs of group 2 with helicity (+,+,-) and (+,-,+) 
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Numerators of the graphs of group 4 with helicity (+,+,-) and (+,-,+) 
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Numerators of the graphs of group 6 with hehcity (+,+,-) and (+,-,+). 



Numerators of the graphs of group 7 with hehcity (+,+,-) and (+,-,+). 
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Numerators of the graphs of group 8 with hehcity (+,+,-) and (+,-,+) 



Numerators of the graphs of group 9 with hehcity (+,+,-) and (+,-,+) 



Numerators of the graphs of group 10 with hehcity (+,+,-) and (+,-,+). 



Basic input parameters. 
Charmonium parameters. 

Baryon wavefunction parameters derived in ref. with the constituent strange 
quark mass = 350 MeV, the octet baryon decay constant /^g = 6.64 x 10~^ 
GeV^ and the transverse size parameter a^g = 0.75 GeV~^ at the reference scale /iq. 
The same for the decuplet baryons are /^^g = 0.0143 GeV^ and a^io = 0.80 GeV^^ 
Note that in ref. [|ll[] of all the decuplet baryons, only the parameters of A were 
given. 
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